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ABSTRACT 

We have conducted an archival XMM-Newton study of the bright X-ray point sources 
in 32 nearby galaxies. From our list of approximately 100 point sources, we attempt 
to determine if there is a low-state counterpart to the Ultraluminous X-ray (ULX) 
population. Indeed. 16 sources in our sample match the criteria we set for a low-state 
ULX, namely, Lx > 10 38 ergs _1 and a spectrum best fit with an absorbed power law. 
Further, we find evidence for 26 high-state ULXs which are best fit by a combined 
blackbody and a power law. As in Galactic black hole systems, the spectral indices, T, 
of the low-state objects, as well as the luminosities, tend to be lower than those of the 
high-state objects. The observed range of blackbody temperatures is 0.1-1 keV with the 
most luminous systems tending toward the lowest temperatures. We also find a class 
of object whose properties (luminosity, blackbody temperature, and power law slopes) 
are very similar to those of galactic stellar mass black holes. In addition, we find a 
subset of these objects that can be best fit by a Comptonized spectrum similar to that 
used for Galactic black holes in the “very high” state, when they are radiating near the 
Eddington limit. 

Subject headings: galaxies: general — surveys — X-rays:binaries — accretion, accretion 
discs 
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1. Introduction 

Through X-ray observations of nearby galaxies, a class of Ultraluminous X-ray (ULX) sources 
has emerged. These are pointlike, non-nuclear sources with bolometric luminosities in excess of 
the Eddington limit for a 20 M© black hole, or L > 2.8 x 10 39 ergs -1 . The true nature of these 
sources is unclear, and this class more likely includes several different types of objects. Though a 
number of these sources are located within a few parsecs of their host galaxy’s dynamical center, 
they do not exhibit many of the characteristics of active galactic nuclei (AGN). Because the ratio 
of X-ray to optical flux is a factor of 10 greater than that of AGN (Anderson et al. 2003; Stocke et 
al. 1983), these objects are fairly easy to recognize in X-ray imaging data. 

Assuming that the Eddington limit is obeyed by black hole accretion, the existence of such 
luminous non- AGN sources presents a puzzle. Several models have been proposed to account for 
the high luminosities of the ULXs. Among these are relativistic and non-relativistic beaming from 
stellar-mass black hole systems (Kording et al. 2002) and accretion of matter into intermediate 
mass black holes (IMBH). In several systems (NGC 1313 X-2, M81 X-9, etc.), detection of emis- 
sion nebulae surrounding the ULX supports isotropic emission from the central source (Pakull 
Sz Mirioni 2003), which cannot be described through relativistic beaming. Further, a number of 
ULX (NGC1313 X-l, etc.) X-ray spectra are best fit with combined multi-component blackbody 
(MOD) and power law fits, similar to Galactic black holes in their high-state. Recently, Miller 
et al. (2003) find that many spectral fits of ULX require cool accretion disk temperatures of ap- 
proximately 100 eV. The theoretical relationship between black hole mass and disk temperature (T 
oc M” 1 / 4 ) has been observed to hold true for stellar mass (typically around 1 keV) and supermas- 
sive (around 10-100 eV) black holes (Makishima et al. 2000). Thus, the cool accretion disk ULXs 
would correspond to a population of high-state IMBHs with masses of « 16 — 10 4 M©. 

If some ULXs do indeed represent a class of high-state IMBHs, similar to the high-state 
stellar mass black holes in our galaxy, we might also expect to see the low-state objects from this 
same population. In Galactic black hole systems, the low-state is generally characterized by lower 
luminosity, with L< 0.1 L Edd (Done & Gierlinski 2003), and a power law spectrum, typically with 
index T « 2.0 (McClintock & Remillard 2004). In this study we seek to find these low-state sources, 
classify the properties of both high-state and low-state ULXs, and examine whether these data are 
consistent or inconsistent with the IMBH hypothesis. 

We present -the results of a detailed analysis of ULXs in nearby galaxies observed with the 
European Space Agency’s XMM-Newton observatory. Only XMM-Newton provides’ the count rates 
and bandpass necessary to distinguish different spectral models for the ULX and accurately deter- 
mine both the temperature of the thermal component expected for high-state objects and whether 
this component is required in the spectral modeling of these objects. 

In Section 2, we detail the observations examined from the XMM-Newton archives and explain 
the data analysis for the individual point sources. In Sections 3 and 4, we discuss the spectral 
fitting technique as well as simulations we conducted to determine their validity. We discuss the 
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implications of our results in Section 5. 


2. Observations and Data Reduction 

The data used in this investigation were drawn from the XMM-Newton public data archive. 
Assuming that low-state ULXs exist in the luminosity range of 10 38 ” 39 ergs -1 , we conducted simu- 
lations to determine the optimum criteria for observations capable of resolving point sources of this 
luminosity. This luminosity range was chosen on the assumption that an approximately 100 Mq 
black hole would radiate at as 5% of the L Edd m the low-state (Done & Gierlinski 2003). Our 
simulations sought to determine the number of photons required to distinguish between spectral 
fits using a power law model, a bremsstrahlung model, and a combined blackbody and power law 
model at an Lx ~ 10 38 ergs -1 , in addition to limiting the possibility of source confusion. We 
adopt these models since they qualitatively correspond to the spectra of low-state black hole X- 
ray binaries, neutron star X-ray binaries, and high-state black hole X-ray binaries, respectively. 
In order to distinguish between the different spectral fits for objects with Lx ~ 10 38 ergs -1 , our 
selection includes galaxies that were observed for at least 10 ks (with the exception of the bright 
ULX in NGC 5408, which had enough photons for analysis despite the low exposure time) with 
XMM-Newton and that are no more distant than 8 Mpc. This yields a minimum of 400 counts, 
objects with Lx > 2 x 10 38 ergs -1 . 

- Our sample of galaxies is selective in that it represents objects of interest in the X-ray band. We 
include details on these host galaxies in Table 1. In most cases, the original choice to observe these 
galaxies (and hence their inclusion in our sample) was unlikely to be influenced by any previous 
knowledge of their ULX population. Hence, this sample should not be biased in terms of the ULX 
population. Our galaxies include only spirals and irregulars. Figure 1 displays the distribution of 
galaxy type. 

We reduced the data using the XMM-Newton Science Analysis System (SAS) version 6.0.0. 
Since the processed pipeline products (PPS) were created with earlier versions of SAS, the ob- 
servation data files (ODF) were used to produce calibrated photon event files for the EPIC-MOS 
and PN cameras using the commands emchain and epchain. Following this, the events tables were 
filtered using the standard criteria outlined in the XMM ABC Guide. For the MOS data (both 
MOSl and MOS2 cameras), good events constitute those with a pulse height in the range of 0.2 
to 12 keV and event patterns that are characterized as 0-12 (single, double, triple, and quadruple 
pixel events). For the PN camera, only patterns of 0-4 (single and double pixel events) are kept, 
with the energy range for the pulse height set between 0.2 and 15 keV. Bad pixels and events too 
close to the edges of the CCD chips were rejected using the stringent selection expression “FLAG 
== 0”. Time filtering was applied as needed by editing the light curve produced in xmmselect for 
the entire observation. Flare events (distinguished by their high count rate) detected in all three 
cameras, were cut using the tabgtigen task as outlined in the ABC Guide. Such filtering was only 
done as needed. 
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Before extracting spectra of the brightest sources, contour maps of the X-ray observation 
were overlaid on Digital Sky Survey (DSS) images. This ensured that bright foreground stars and 
background AGN were easily distinguished, and thereby not included in the spectral fitting. Also, 
we checked the XMM-Newton positions with NED and SIMBAD to determine if they coincide with 
any known background galaxies or QSOs. A list of these bright fore-ground or background sources 
is included in Table 7. 


3. Spectral Fitting 

Spectra for the bright point sources were extracted using the SAS task especget With this task 
we created spectra (for both the source and background), response matrices, and ancillary response 
files for all three EPIC cameras, when possible. The typical extraction radius was 20 arcseconds, 
but depending on both the size and proximity of a source to another source, the extraction radius 
ranged from 9-87 arcseconds. Background spectra were extracted either in an annulus centered 
on the source, or in a circle of appropriate size away from the source, depending on the proximity 
of the central source to other X-ray sources. Once the spectra were obtained, they were rebinned to 
require at least 20 counts per bin, using the command grppha in LHEASOFT. The list of sources, 
with position and count information, is included in Table 6. 

The extracted spectra were fit with standard models in XSPEC vl 1.3.1. For each source, we fit 
the PN and MOS spectra simultaneously in the 0.3-10 keV range. We allowed a free normalization 
constant to account for the differences in flux calibration between the three cameras (similar to 
Jenkins et al. (2004)). Each source was first fit with an absorbed single component model. In all 
cases we used the standard absorption model wabs , leaving the column density as a free parameter. 
Results of the single-component fits are seen in Table 2. We include in this table only the best-fit 
parameters for those sources best described by a single-component model. The flux values quoted 
represent the unabsorbed flux in the PN spectra, in the 0.3-10 keV band. All errors quoted, here 
and subsequently, correspond to the 90% confidence level for one degree of freedom (A% 2 = 2.71). 
The luminosities were calculated from the unabsorbed flux using the distances quoted in Table 1. 
Both flux and luminosity correspond to those of the best fit model (power law or bremsstrahlung). 
It should be noted that since our selection criteria was based on a count rate cutoff, due to the 
variety of spectral forms, the inferred luminosity cutoff will not be uniform. 

For a number of sources, the single-component models did not adequately describe their spec- 
tra. For these sources we employed an absorbed two-component blackbody and power law model. 
In Table 3 we present the results of the sources best fit by this two-component fit. We include 
the improvement in x 2 of the two-component fit over the simple power law. We include the power 
law best fits to these sources in the appendix. Table 4 includes those sources where both the 
single-component and two-component models were indistinguishable. 

Further, for those sources we classify as ULXs (see section 5 for our criteria), we computed 
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bolometric luminosities. We used the exponentially cutoff power law spectrum of Magdziarz & 
Zdziarski (1995), model pexrav in XSPEC, with a cutoff energy of lOkeV and the reflection pa- 
rameter set to zero. This model was used in place of the power law component. We computed an 
unabsorbed flux in the 0.1 - 100 keV range through use of the dummyresp command (which extends 
the model beyond the observation’s energy range). The luminosity was then computed using the 
distances listed in Table 1. We quote these values as L upper in Table 5. We note that these values 
represent an upper limit on the bolometric luminosity for steep power law (T > 2) objects, since we 
would expect the power law component to cutoff at some low energy. However, for flat spectrum 
(T < 2) sources L boi is a lower limit. 

For our ULX sources modeled by a combined blackbody and exponentially cutoff power law, 
we estimate a more accurate L(, 0 ;, calculated from the flux in the range of 2x kT- 100 keV where 
kT is the blackbody temperature obtained from the model. In galactic X-ray binary systems, the 
power law component of the X-ray spectrum is believed to be from Comptonization in a corona. 
The photons, supplying this energy originate from the blackbody continuum emanating from the 
accretion disk. Thus, a natural cutoff for this power law component occurs at the peak emission of 
the blackbody (which is approximately 3x kT). These estimated blackbody values are within 95% 
of the full integrated blackbody flux and are therefore a good approximation to the data. 

We note that our bolometric luminosities, on average, are a factor of 1.08 greater than the X- 
ray luminosities in the 0.3 - 10 keV band for the objects best fit by a combined blackbody and power 
law. Thus, to good approximation, the X-ray luminosity is the bolometric luminosity. However, 
for the objects best fit by a simple power law, the average bolometric luminosity is roughly a factor 
of 7 greater than the X-ray luminosity in our band. This average is dominated by the steep power 
law objects, in particular Holmberg II XMM1 (T = 3.09). Excluding this object, we get an average 
bolometric luminosity that is 2.8 times the X-ray flux and more indicative of the general properties 
of these power law-fit objects. 

In this large sample of point sources, we came across a number of objects whose spectra were 
not well fit by the standard models employed. We briefly describe these sources in the appendix. 


4. Spectral Simulations 

In order to determine whether the blackbody component is statistically significant for the 
sources fit with a two-component model, we simulated spectra of some of the brightest sources. 
Table 8 shows the results of these simulations. We simulated spectra using the XSPEC command 
fakeit for 10 bright (> 5000 counts) sources that were best fit by a combined blackbody and a 
power law model. The fakeit command uses Monte Carlo simulations in order to create a simulated 
spectrum from the original dataset. We simulated spectra with both 10% and 5% of the actual 
exposure time. These spectra were then fit with the three standard spectral fits employed in this 
study. 
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Comparing chi-square values for the simulated spectral fits, the combined blackbody and power 
law fit is the best fit for all simulated spectra at 10% the actual exposure time, corresponding to a 
source that is ten times weaker. At 5% the exposure time, 70% of the simulated spectra are well-fit 
by the combined model. These results suggest that indeed the blackbody component is statistically 
significant and distinct from a pure absorbed power law model. Thus, we are confident that the 
best fit for the bright combined blackbody and power law spectral fits is a valid description of the 
data. 


5. Discussion 

We have determined best-fit spectral parameters of the bright X-ray sources in 32 nearby 
galaxies. In choosing three “standard” models for our study, we hoped to accurately separate high 
and low state ULXs from other bright X-ray sources. We specifically chose to fit the data with the 
bremsstrahlung model in order to identify neutron star X-ray binaries within our sample. 

We cross-referenced the X-ray positions of our sources with both NED and SIMB AD in order to 
identify known supernovae, galaxies, and stars. In addition, we examined the DSS optical images to 
place the position of our sources within their respective galaxies. Such analysis aimed to minimize 
contamination of our sample of ULXs with bright background and foreground sources. 

Further, we examined XMM-Newton’s Optical Monitor data in the visual bands (U, B, V). The 
PPS contain point source detection files for the OM data. We overlaid these point source detections 
with X-ray contour maps in order to determine the brightest possible optical count rates for the 
X-ray sources, which were then converted into fluxes using the OM calibration documentation. In 
Figure 2, we plot the distribution of the logarithm of the X-ray to optical flux for the brightest 
possible optical counterpart inside the XMM-Newton error circle. Only 13 of the 32 host galaxies 
had visible OM data during the observations. Of these 13 galaxies, 40 sources were in the range of 
the OM data and only 14 were coincident with an optical point source. Therefore, the majority of 
our sources have X-ray /optical flux ratios that are larger than those displayed. We estimate the 
point source detection limit of the OM U filter as approximately 1.24 x 10” 14 erg cm” 2 s” 1 . For an 
unabsorbed X-ray flux of 1. Ox 10” 12 erg cm” 2 s” 1 , this corresponds to log(/ x // op f) — 1.9. Therefore, 
the average value for our sources should fall around 2 or greater. The average distribution for QSOs 
and AGN centers around 0 and 0.8 for BL Lacs (Anderson et al. 2003). Thus, our objects have 
ratios of L x /L opt atleast 10 times higher than those of AGN and 100 times greater than stars. 

Recently, Gutierrez & Lopez- Corredoira (2005) identify six ULXs from the catalog of Colbert 
& Ptak (2002) as QSOs. They hypothesize that a large number of ULXs may in fact be quasars at 
higher redshift than their supposed host galaxy. However, unlike the objects studied in Gutierrez Sz 
Lopez-Corredoira (2005), our ULX sources are all spatially coincident with the optical host galaxy. 
In addition, a majority of our ULXs are not in the proximity of a cataloged optical point source. 
The X-ray/optical flux ratios of our sources are much larger, on average, than might be expected for 
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a QSO. It is also worth noting that while some cataloged ULXs may be QSOs, optical identifications 
have been made associating other ULXs with a type B supergiant companion (Kuntz et al. 2005; 
Liu et al. 2004). 

5.1. Classification Criteria 

The spectral fits indicate that we can indeed distinguish a class of low-state ULXs from the 
high-state objects. As a constraint on our ULX classification, we require our X-ray sources to 
clearly coincide with the optical extent of the host galaxy (as determined from the DSS images) . 
Of the sources in Table 2, 16 are “low-state” objects, having unabsorbed luminosities > 10 38 ergs -1 
and spectra that are best fit by power law models. Further, 26 sources have unabsorbed Lx > 
3 x 10 39 ergs -1 and spectra that are best fit by combined blackbody and power law models. These 
are “high-state” objects. In addition to these high and low state ULXs, we find a large number 
of sources best fit by a combined blackbody and power law model but below our threshold of 
Lx = 3 x 10 39 ergs -1 . Many of these sources may be accreting stellar mass black holes. Some of 
these sources were found away from the optical extent of the targeted galaxy (from our analysis of 
the DSS images), and may represent background AGN. 


5.2. Low-State ULX 

For Galactic black hole X-ray binaries, spectral indices of low-state objects are typically lower 
than those of high-state objects (McClintock & Remillard 2004). In Figure 3, we plot the distri- 
bution of the spectral index for both high-state and low-state objects. As in the Galactic sources, 
it is clearly shown that the spectral indices of the high-state objects are indeed larger. Of further 
interest, the distribution of spectral index for low-state objects looks remarkably similar to the 
distribution of spectral index for moderate luminosity quasars, many of which are thought to be 
the analogs of low-state black holes (Porquet et al. 2004). This supports the classification of these 
objects as accreting black holes. 

The low hard X-ray state of X-ray binaries is associated with a low accretion rate from the 
companion object. Therefore, on average, we expect the luminosities of the low-state objects to 
be lower than the high-state objects. Figure 4 displays luminosity as a function of the spectral 
index. On average, the highest luminosity low-state objects have luminosities lower than those of 
the high-state objects. We find mean values of T = 2.46, with a root mean square (rms) deviation 
of S= 0.12, and Lx = 1-4 x 10 4 °ergs -1 , log(S) = 1.6, for the high-state objects. This calculation 
excludes the 3 objects with spectral indices greater than 3.5. For the low-state objects, we find 
mean values of T = 2.09, with a rms deviation of S= 0.10, and Lx = 2.2 x 10 39 ergs -1 , log(S) = 2.1. 

The lower Lx values of the low-state objects imply that they may indeed be accreting at a 
lower rate than the high-state objects. This can further be seen in the bolometric luminosities 
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listed in Table 5. If these objects are accreting at a rate similar to galactic low-state black holes 
(0.1 x L£dd) (Done & Gierlinski 2003), we can estimate their masses as 

M _ Lboi 

Mq 0.1 X L E dd 

with L E dd as the Eddington luminosity for a 1 M 0 object (1.3 x 10 38 ergs -1 ). Our mass estimations, 
based upon our limits to the bolometric luminosities, yield masses of 20 - 1524M 0 (see Table 5), 
precisely what we might expect for a population of IMBHs. 


5.3. High-State ULX 

If the high-state ULXs represent a class of intermediate mass black hole systems, their X-ray 
spectra should be best fit by a combined blackbody and power law model. Scaling for the mass of 
the black hole, we would expect a relationship of T oc M -1 / 4 between black hole mass and blackbody 
temperature (Makishima et al. 2000). This would indicate a thermal component of ~ 100 eV. A 
few objects have been reported to display this property (Miller et al. 2003; Roberts Sz Warwick 
2000). In Figure 5, we graph the distribution of the thermal component for our classified high-state 
objects. 

We find that there are two peaks in the distribution among the thermal component, one at 
approximately 100 eV and another centered close to 1 keV. This could indicate two different classes 
among the high-state objects. It is possible that those objects with blackbody components near 100 
eV are indeed high-state intermediate mass black holes. The second peak, centered around 1 keV, 
has thermal components reminiscent of the Galactic black hole systems in our own galaxy. These 
systems may be stellar-mass black holes accreting matter at an increased rate. If this were the case, 
we would expect the luminosities of the sources exhibiting a higher blackbody temperature to be 
lower. In the second graph of Figure 5, we plot the relationship between blackbody temperature 
and Lx in the 0.3 - 10 keV band. Once again, two classes of ULXs are seen. The most luminous 
objects are those with low blackbody temperatures. On average, the less luminous sources exhibit 
higher blackbody temperatures. 

The second, low-luminosity, class of ULX is clearly distinguishable in both plots of Figure 5. 
We found that, with the exceptions of NGC 253 XMM1, M81 XMM1, and NGC 5204 XMM1, the 
spectra of these objects could be well-described by an absorbed Comptonization ( compST) model 
(Sunyaev & Titarchuk 1980) used to fit galactic black holes in the “very high” state when they 
are radiating at the Eddington limit. This model simulates Compton scattering of cool photons 
on the hot electrons of a completely ionized plasma. We present the best-fit parameters for the 
Comptonization model in Table 10. 

This “very high” state has been observed (Miyamoto et al. 1991) in a few Galactic black holes. 
McClintock & Remillard (2004) use the alternative nomenclature of the steep power law state, a 
state which is characterized by T > 2.4 and a luminosity which may or may not be greater than 
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the luminosity in the high-state. Yet another rubric for the very high state emerged in Kubota et 
al. (2001) and Kubota & Makishima (2004), where they identify this as the “anomalous” state, a 
state whose spectrum can be well fit by an inverse Compton scattering model. Regardless of the 
name, our best-fit Comptonization sources likely fit into this category. The luminosities of these 
sources suggest that they are stellar mass black hole systems in this anomalous/ very high state. 

As with the low-state, we include mass estimates for our high-state objects in Table 5. We 
assume that the high-state objects are radiating at L Edd- We find masses of 1.6 — 38 M© for the 
sources well fit by the Comptonization model. The other high-state ULXs range from 16.5 - 
1354 M©, analogous to the low-state ULX masses computed. 


5.4. Temperature Gap 

In addition to the existence of ULXs with low blackbody temperatures, the temperature distri- 
bution (Figure 5, left panel) displays a “gap” which is of particular interest — there is a complete 
absence of objects with temperatures in the range 0.26 keV to 0.50 keV. It is tempting to take this 
as evidence for a gap in the mass distribution of these accreting black holes. Since, for a given 
luminosity, we expect the temperature to vary as T oc L 1 / 4 M” 1 / 2 , this factor of two gap in the 
temperature distribution translates into a factor of four gap in the black hole mass distribution. 

If this result is borne out by further study, it provides an important clue to the origin and 
evolution of intermediate mass black holes. One popular idea is that intermediate mass black holes 
formed from the collapse of massive Population III stars (Madau & Rees 2001). Models suggest 
that Pop III stars with zero age main sequence (ZAMS) masses in the range 25-140M© and above 
260M© collapse to produce black holes (Heger & Woosley 2002) whereas in the range of ZAMS 
masses 140-260M©, pair-instability supernovae lead to the complete disruption of the stars (i.e., 
no remnant black hole remains). Hence, this model for IMBH formation predicts a gap in the 
IMBH initial mass function in the range of approximately 60-200M© (although this is uncertain 
on the low end due to the effect of the pulsational pair- instability on the pre-collapse core). One 
possibility is that the gap in our observed temperature distribution (and hence the inferred gap 
in the mass function) is due to this effect of the pair instability supernovae in Pop III stars. This 
would require that the current IMBH mass function is approximately the same as the initial IMBH 
mass function. In other words, it requires that most IMBHs (especially those just below the gap) 
have not grown significantly due to accretion since their formation and, hence, that the ULX phase 
represents a short fraction of the life-time of an IMBH (/ << isalAm where t sa \ ~ 45eo.i Myr is the 
e- folding timescale for Eddington limited black hole growth with radiative efficiency e — O.leo.i)- 

An alternative interpretation of the inferred mass gap is to suppose that two fundamentally 
different modes of formation lead to a strong bi-modality in the final black hole mass function. 
Black hole masses below the gap can be readily understood through normal stellar processes. A 
separate and distinct population of significantly more massive black holes may result from dynamical 
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processes in the core of dense globular clusters (Miller & Hamilton 2002; Gultekin, Miller, Sz 
Hamilton 2004). 


5.5. Galactic HMXBs 

Supposing that the Galaxy’s bright X-ray population is representative of low-redshift galaxies, 
we expected to find a number of sources similar to Galactic X-ray binaries in our sample. As 
previously stated, we set a luminosity cutoff of x 10 39 ergs~ 1 (0.3 — lOkeV band) in order to 
distinguish between galactic HMXBs and high-state ULXs. In our sample, we find approximately 
24 sources with luminosities below our high-state ULX cutoff, X-ray positions within the optical 
extent of their host galaxy, and no obvious optical counterpart. The unabsorbed luminosities for 
these sources range from 0.4 - 2.5 x 10 39 ergs“ 1 (0.3 — lOkeV band). Two of these sources were 
transients. Of the four host galaxies with multiple observations examined, two of these galaxies 
contained solely ULX sources in our luminosity regime (Holmberg II and NGC 5204). Each of the 
remaining two (NGC 253 and NGC 4258) had a transient source best fit by a combined blackbody 
and a power law. 

This suggests an interesting diagnostic in terms of distinguishing our ULX sources from a 
normal HMXB population. In our own galaxy, most HMXBs vary on timescales of days or less 
and most of the black holes in the Milky Way are transients. The figures in a recent paper 
of Kalogera et al. (2004), determined through detailed mass- transfer calculations, indicate that 
transient behavior should not be expected from a population of IMBHs (M. Coleman Miller 2005, 
private communication). Thus, on average, our ULX sources should remain X-ray bright in multiple 
observations. Through a literature search, we found that 37/42 ULX sources are distinguishable 
in ROSAT observations and thus are luminous for greater than 10 years and therefore are not 
transients. Examination of the long term light curves show that most sources vary by less than a 
factor of 3 over the timescale from ROSAT to XMM. 

As a possible further diagnostic, we constructed a color-color diagram for our ULX sources. 
We adopted the colors of Done & Gierlinski (2003) in order to compare our sample with their 
sample of Galactic X-ray sources. Thus, our colors were constructed from unabsorbed model fluxes 
in four energy bands: 3-4, 4-6.4, 6. 4-9. 7, and 9.7-16 keV. The XSPEC command dummyresp was 
used to calculate a flux based on the model for the 10-16 keV range. We plot colors for a pure 
unabsorbed power law (from T = 1.5 — 3.0) and an unabsorbed MCD model ( diskbb in XSPEC with 
kTi n = 5.0 — 0.2 eV) for comparison. Comparing our Figure 6 with Figure 8 of Done & Gierlinski 
(2003), we find that our ULX sources largely lie along the same regions as their black hole sources. 
A few ULX sources, however, lie in the region occupied by atoll and Z-sources in the plot of Done 
& Gierlinski (2003). These sources were those best fit by a Comptonization model. 
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5.6. Galaxy Sample 

It has been shown that the ULX population is proportional to the host galaxy’s star formation 
rate (SFR) (Ranalli et al. 2002; Grimm et al. 2003). The far-infrared luminosity of a galaxy is used 
as an indicator of the SFR. In order to compare the ULX population of a galaxy with the SFR we 
followed a similar approach to Swartz et al. (2004). We calculate the FIR flux from observations 
taken by the Infrared Astronomical Satellite. As in Swartz et al. (2004), the flux between 42.4 
and 122.5 fim is approximated as 1.26xlO -u (2.58S6o+Sioo) ergcm -2 s -1 . The values of the flux 
at 60 ^m (Seo) and 100 nm (Sioo) were obtained from either Ho et al. (1997) or NED. Luminosities 
were calculated using the distances quoted in Table 1. We list these values in addition to the 
number of ULXs observed in individual galaxies in Table 11. The number of ULXs includes both 
the objects we classify as high and low state ULX as well as those sources resolved by Chandra. 

In Figure 7, we show two plots relating the number of ULXs to Ijfir- It has been suggested 
by Grimm et al. (2003) that the luminosity function in the X-ray regime from HMXBs is related 
to SFR. In our first plot, we find that the galaxies with the highest Lfir seem to have fewer ULXs 
than may be expected from the luminosity functions of Grimm et al. (2003). Thus, in a direct 
comparison, our results do not agree with their predictions. 

The second plot displays the average number of ULXs /galaxy, binned according to luminosity. 
This plot is extremely similar to Fig. 15 of Swartz et al. (2004) for spiral galaxies. Thus, once 
again, it seems that the connection between SFR and the ULX population in spirals is supported. 
For irregular galaxies, however, there seems to be more of a spread in the distribution. This could 
be the result of poor sampling — most of the bins contain only one galaxy. Another possibility is 
that there is no direct correlation in irregular galaxies or that the overall star formation in these 
galaxies is less ordered or clumpier. If the latter is the case, the overall SFR of the galaxy is only 
an average over a wide range of values. We shall address this issue again in the next paper in this 
series (L.M. Winter et al., in preparation) where we discuss the local environment of the ULX in 
our sample. 

In Figure 8 we plot the distribution of column densities among the ULX. We subtracted the 
Galactic column density towards the galaxy (obtained from the nH FTOOL and listed in Table 1) 
from the values obtained through spectral fits. We note that on average the ULXs have large 
column densities. The typical Galactic column density along a line of sight is « 4 x 10 20 cm -2 . 
If the ULX is located on the opposite side of its host galaxy, we might expect maximum column 
densities of « 1.2 x 10 21 cm -2 . However, most of our sources have column densities well above 
this value. This may imply that the local environment of the ULXs contains an extra source of 
absorption. 

In order to better understand the relationship between SFR and the ULX population, it is 
necessary to extend ULX studies to other wavelengths. In particular, it becomes important to 
analyze UV and IR images close to the ULX. 
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6. Conclusion 

We have found from our XMM survey of the ULX population in nearby galaxies that there 
exists a population of objects whose X-ray spectral properties closely match the low-state spectra of 
Galactic black holes, but whose luminosities lie in the range of L ^ - 2 x 10 38 -lx 10 40 ergs -1 . In 
the Milky Way, black holes with these spectral properties radiate at only ^ 0.05 of the Eddington 
limit. If this is also true for this population, it indirectly implies that these objects have a mass 
greater than ^ 30 M 0 ranging up to 1500 M 0 and thus should be IMBHs. The existence of such 
objects was “predicted” on the basis that the ULXs previously studied shared the X-ray spectral 
characteristics of high-state Galactic black holes; namely, an X-ray spectrum best fit by a combined 
blackbody and a power law (Miller et al. 2003), but with much higher luminosities. If these objects 
are high-state IMBHs, the corresponding low-state objects should also exist. 

Our survey has also uncovered a large population of objects whose X-ray spectra are well mod- 
eled by the canonical description of Galactic black holes in the high-state, a black hole with a steep 
power law, but whose bolometric luminosities exceed 2 x 10 39 ergs -1 , ranging up to 10 41 * 5 ergs -1 
and whose blackbody temperatures are less than 0.3 keV. If these objects are radiating at « 1/2 
the Eddington limit like their Milky Way counterparts their implied masses are from 30 — 3000 M 0 , 
a range very similar to that implied by the low-state objects. Using the M -1 / 4 scaling of mass 
to temperature, the observed spectral temperatures give masses of 500 - 10 4 M 0 a considerably 
larger value. In general agreement with the expectations of the IMBH hypothesis, the objects with 
high-state spectra are more luminous than those with low-state spectra. 

We conclude, from an X-ray spectral and luminosity point of view, that our data are consistent 
with many of these objects having the properties expected of an IMBH population. However, we 
also find two other populations of objects, those whose blackbody temperature and luminosity 
correspond to that of stellar mass black holes with kT ~ 1 keV and log Lx less than 2 x 10 39 ergs -1 
and a small population of objects whose X-ray spectra and luminosities are consistent with that 
of stellar mass black holes in the very high state. Thus, ULX selected purely on the basis of 
0.3- lOkeV X-ray luminosities are a composite class with « 1/4 being “normal” stellar mass black 
holes and the rest being consistent with a population of IMBHs. 

The existence of a substantial population of ULXs in nearby dwarf and other low star formation 
rate galaxies argues that (in agreement with Ptak & Colbert (2004); Swartz et al. (2004)) there is 
more than one source term for the origin of ULXs, with at least some of them not being associated 
with recent star formation, at least statistically. We note that these results have required the high 
signal to noise of XMM in order to discern the spectrum of these objects. Many of these objects 
have also been observed by Chandra and their spectra have been well-fitted by simple power laws. 

In a follow-up paper we will discuss the environments of these objects as revealed by XMM 
OM UV imaging and the implications this has for the origin of ULXs. 
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A. Appendix material 

The following sources were not best fit by the standard models employed in this study: 


A.l. NGC300 XMM4 

This source was classified as a super-soft X-ray source by Kong & DiStefano (2003). We find 
that the standard single-component absorbed blackbody model is a much better model for this 
spectrum. In fact, the power law, bremsstrahlung, and combined models do not fit the data within 
the 90% confidence range. Fitting an absorbed blackbody, we find the best fit corresponds to the 
following parameters: n# = 1.38i°;55 x 10 21 cm -2 , kT = 0.0591q;qo 5, and x 2 /dof = 74.5/45. This 
fit yields an unabsorbed flux of 3.3 x 10~ 13 ergcm -2 s -1 . 


A. 2. NGC4631 XMM4 

The spectrum of this source clearly identifies it as a super-soft X-ray source. As with NGC300 
XMM4, the standard models employed in this study did not adequately match the data. The best 
fitting model corresponds again to an absorbed blackbody. The corresponding parameters are as 
follows: tih = 6.21J; 26 x 10 21 cm -2 , kT = 0.07+o;°{, and x 2 /dof = 142.3/74. This fit yields an 
unabsorbed flux of 9.5 x 10 -12 erg cm -2 s -1 . The position of this source shows it to be coincident 
with a globular cluster associated with that galaxy. This source was identified as a bulge X-ray 
source, possibly powered by accretion, in a ROSAT study of NGC4631 (Vogler Sz Pietsch 1996). 


A. 3. NGC4631 XMM5 

The spectrum of this source was best fit with an absorbed power law + an absorbed vapec 
model. This indicates the prescence of hot gas but we have no further explanation. 


A. 4. NGC4945 XMM5 

The spectrum of this source was not adequately fit with any of the standard models used in 
this investigation. The spectrum exhibits a prominent Fe K line in the PN spectrum that is well fit 
by a gaussian ( zgauss ) at 6.4 keV. We find that the entire spectrum is best fit with a partial covering 
fraction absorption model ( pcfabs ) in combination with the normal absorption, a power law, and 
a gaussian. The best fit parameters yield: absorption column density, n# = 1.79 x 10 21 cm -2 , 
partial covering absorption, n# = 18.4 x 10 21 cm -2 , partial covering fraction = 0.82, T = 1.6, and 
X 2 /do f — 61.8/57. The source is clearly located within the optical galaxy, and is thus unlikely to 
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be a background AGN, but we have no other explanation. 


A. 5. M51 XMM5 

The spectrum and luminosity (Lx ~ 1.9 x 10 42 ergs _1 ) of this source suggests that it is an 
AGN. The location of the source, from the Digital Sky Survey, places it within the dwarf companion 
of M51 making a value of the optical flux hard to constrain. The best fit to this source was an 
absorbed blackbody + power law and the spectral parameters are listed in Table 4. 


A. 6. M83 XMM2 

Like NGC4945 XMM5, this source was best fit by a partial absorption model. However, this 
source showed no evidence of an Fe K line. We fit this source’s spectra using a partial covering 
fraction absorption model in combination with the normal absorption model and a power law. 
The best fit parameters yield: absorption column density, n# = 2.1 x 10 21 cm -2 , partial covering 
absorption, nn = 43.5 x 10 21 cm -2 , partial covering fraction = 0.86, T = 2.95, and y^/dof = 
83.5/84. The unabsorbed flux in the range of 0.3-10 keV equals 1.37 x 10~ 12 erg cm -2 s -1 . 


A. 7. Inverse Compton Scattering Sources 

Table 7 includes the parameters for the "ULX” sources best fit by the compST model. A 
discussion of these sources and interpretation of the data is included in section 5.2. 
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Table 1. XMM-Newton Galaxy Observations 


Galaxy 

Type a 

n H b 

distance 0 

ref 

obs id d 

duration (s) 

comments 

NGC247 

SAB(s)d 

1.54 

3.09 


0110990301 

14536 


NGC253 

SAB(s)c;HII 

1.40 

3.73 


0110900101, 0152020101 

30711, 110591 

Starburst 

NGC300 

' SA(s)d 

3.11 

2.56 


0112800101 

43967 

- 

NGC625 

SB(s)m? sp; HII 

2.15 

2.62 


0085100101 

26288 

- 

NGC1313 

SB(s)d; HII 

4.0 

4.17 


0106860101 

41310 

- 

IC0342 

SAB(rs)cd; HII 

30.3 

3.9 

1 

0093640901 

11217 

- 

NGC1569 

IBm 

21.7 

1.6 

1 

0112290801 

15582 

Starburst 

NGC1705 

SA0- pec; HII 

3.9 

5.1 

2 

0148650101 

58926 

Starburst 

MRK 71 

BCD; HII 

3.9 

3.4 

3 

0141150201 

45919 

galaxy pair 

NGC2403 

SAB(s)cd; HII 

4.15 

3.56 


0150651201 

11415 

- 

Holmberg II 

Im 

3.42 

2.70 


0112520701, 0112520901 

13528, 6860 

- 

Holmberg I 

IAB(s)m 

3.49 

3.6 

4 

0026340101 

26280 

- 

M81 

SA(s)ab;LINER 

4.12 

3.6 

4 

0111800101 

127913 

Hoi IX also in field of view 

M82 

10; HII 

4.14 

3.9 

5 

0112290201 

29387 

Starburst 

Holmberg IX 

Im 

4.0 

3.6 

4 

0112521001 

10350 

M81 also in field of view 

Sextans A 

IBm 

3.85 

1.4 

6 

0026340201 

21618 

- 

IC 2574 

SAB(s)m 

2.29 

3.6 

7 

0026340301 

24263 

bursting star-formation 

NGC 4214 

IAB(s)m; HII 

1.49 

2.7 


0035940201 

14744 

- 

NGC 4258 

SAB(s)bc;LINER 

1.2 

7.2 


0059140901, 0110920101 

16146, 21895 

- 

NGC4395 

SA(s)m;LINER 

1.33 

4 


0112521901 

15842 

- 

NGC4449 

IBm; HII 

1.39 

3.08 


0112521701 

15522 

- 

NGC4490 

SB(s)d 

1.78 

7.8 

1 

0112280201 

17754 

interacting with NGC4485 

NGC4631 

SB(s)d 

1.28 

7.5 

1 

0110900201 

53850 

- 

NGC4736 

(R)SA(r)ab;LINER 

1.43 

4.3 

1 

0094360601 

23461 

- 

NGC4945 

SB(s)cd; Sy2 

159 

3.1 


0112310301 

23062 

- 

NGC 5204 

SA(s)m; HII 

1.42 

4.8 

1 

0142770101, 0142770301 

19205, 16387 

- 

M51 

Sc; Sy2 

1.55 

7.2 

1 

0112840201 

20924 

Galaxy pair 

M83 

SAB(s)c;HII 

3.94 

6.2 


0110910201 

30627 

Starburst 

NGC5253 

Im pec;HII 

3.77 

3.2 

1 

0035940301 

47216 

Starburst 

M101 

SAB(rs)cd 

1.17 

7.4 

8 

0104260101 

43019 

- 

NGC5408 

IB(s)m; HII 

5.73 

4.8 

9 

0112290601 

7757 

- 

Circinus 

SA(s)b; Sy2 

57.8 

4 

10 

0111240101 

110496 

- 


a from the NASA/IPAC Extragalactic Database (NED) 

b column density in units of 10 20 cm -2 , obtained from the web version of the nH FTOOL 
C distance in Mpc (if no reference is given, obtained from the distance modulus given in LED A) 
d XMM-Newton observation ids for the data examined in this survey 


References. — (l)Tully 1988; (2) Tosi et al. 2001; (3) Tolstoy et al. 1995; (4) Freedman et al. 1994; (5) Sakai & Madore 1999; (6) Sakai, 
Madore, &: Freedman 1996; (7)Shapley, Fabbiano, & Eskridge 2001; (8) Kelson 1996; (9) Karachentsev et al. 2002; (10) Freeman et al. 1977. 
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Table 2. XMM-Newton best fit: single component spectral fits 




Powerlaw 



Bremsstrahlung 



Source 

n# a 

r 

X 2 /dof 

•n H a 

kT (keV) 

X 2 /dof 

Fx b 

Lx c 

NGC247 XMM2 

1 4 +1-8 

o oq+ 1-02 

Z - zy -0.57 

47.7/54 

0.54±g 11 s 

9 cc+6.90 
z - 00 -1.61 

48.8/54 

0.33 

0.38 

NGC 253 XMM2 (obs 1 ) 

1 6 - * -0 ' 4 


69.1/74 

n k+ 0-2 
uo -0.3 

9 i 9+0.52 
Z - iZ -0.37 

74.7/74 

0.52 

0.87 

NGC300 XMM4 d 

2.5 

9.07 

90.6/45 

0.27 

0.14 

117.6/45 

- 

- 

NGC1313 XMM4 

1 . 86 lgj 

i O+0.07 
i ’°- 0.12 

141.7/149 

1 9 +O .3 
A * Z — 0.3 

6.62l 2 j 8 

140.1/149 

0.33 

0.69 

IC0342 XMM1 

r O + 0.6 

d -°-0.3 

i- 68 ±g;g| 

159.5/185 

4-9±g;| 

I0.5l?;| 

160/185 

3.5 

6.37 

IC0342 XMM2 

23.9t$:S 

1 ftc;"t 0.22 

l.OO _0 20 

9 Q 9 + 0.20 

z - UZ -0.19 

77.5/85 

2 i±|;l 

8.51?;° 

74.9/85 

4.64 

8.44 

IC0342 XMM4 

5.3111 

1 y+ 0.8 

-0.7 

1 C S + 9 ' 2 
lo - 0.2 

64/58 

a 9+6.99 
4 ' Z -0.85 
1 1 + 0.5 
i ' i -0.3 

4.441J;?? 

56.9/58 

0.69 

1.26 

NGC2403 XMM4 

1 OQ+O - 39 

l-oy_0.25 

62.3/71 

4-591?;} 

62.3/71 

0.31 

0.48 

Holmll XMM1 (obs 2) 

3 0 Q +°* 15 
°*u y _o ,12 

266.7/252 

0 . 3 +!;}! 

1 IS 4-0 * 11 
iA, 3 -o.ii 

309.4/252 

3.5 

3.1 

Holm I XMM2 

0.35 e 

9 i 0 +O.I 6 
Z ' id -0.15 

39.2/45 


9 r+0.79 
Z '°-0.55 

51.6/45 

0.10 

0.16 

Holm I XMM3 

0.35 e 

2.051°;}! 
1-97 ±g;?5 

34.4/32 

0.35 e 

9 ao+0.85 
z * u ^ — 0.56 

42.1/32 

0.12 

0.19 

IC2574 XMMl 

i o+0.40 
1 — 0.30 

120.9/103 

0.69i°; 2 ! 

4 1+0-89 

107.5/103 

0.35 

0.47 

IC2574 XMM2 

0 4 “*" 0 ‘ 4 
u -^-0-3 

0 9 + 0.21 
zz -0.09 

45.7/51 

?t? 

1-971?;? 

57/51 

0.22 

0.34 

NGC4214 XMMl 

1 1 +0.52 
li -0.47 

1 Q7+0-26 

A ' 5 ' - 0.21 
9 4Q+0-44 
z *"*»_0.33 
9 09 +O .34 
z *^-0.24 

41.9/38 

0.54t°; 4 ! 

4.861?;?! 

2 61+ 1 - 22 
Z ’ Oi -0.72 

44.5/38 

0.25 

0.22 

NGC4258 XMM2 (obs 2 ) 

6 7 +2 - 6 
D,/ -1.5 

83.6/57 

4.8i°;l 

85.5/57 

0.30 

1.9 

NGC4258 XMM3 

1 4 + 0-69 
i - 4 -0.64 

38.9/37 

o.49i!i! 

9 40+ 1 - 09 
^•4S — o.74 

41.3/37 

0.20 

1.2 


3.8 

1.82 

4/11 

2.7 

7.14 

5/11 

0.077 

0.48 

NGC4258 XMM4 

0 . 68 i°; 2 J 

1 Q 7 + 0.22 
A * y ' -0.19 

41.1/48 

ooeto.il! 

4.07l};| 

45.2/48 

0.39 

2.4 


i n + 0.78 

i-y_o. 6 o 

2 24~ t "°‘ 29 
ZlZ4 -0.24 

77.03/77 

o. 9 lg;S 

2.821J;! 

77.8/77 

0.33 

2.0 

NGC4395 XMM2 

0.33i“j 

9 7 C+O .45 
z *'°-0.33 

38.6/36 


0.891°;}? 

50/36 

0.15 

0.28 

NGC4395 XMM4 

°-3tg; S 

2-08lg;|g 

16/25 

?±t 

9 97 +I.I 6 

Z * Z i — 0.72 

19.9/25 

0.15 

0.28 

NGC4449 XMMl 

6 3 +°- 9 

0.^-0 7 

9 99+0.14 
ZlZZ - 0.12 

103/118 

4-5±g;| 

0 Q7+0.73 
5 ‘ y ‘ -0.58 

1 fi ^+ 0,22 
l.DO_o 21 

114.3/118 

1.2 

1.36 

NGC4449 XMM2 

1 5 +0 - 3 
1 °-0.3 

9 qi +0.16 
Z '° A -0.14 

103.5/112 

0.25t°; 2 

112.1/112 

0.29 

0.33 

NGC4490 XMM4 

10.217.1 

2 0 Q + 0 - 23 
zuy -0.19 

51.6/50 

8.3+};! 

4-751!;!! 

3.08l°;!? 

50.3/50 

0.84 

6.1 

NGC4490 XMM5 

o q + 0.94 

9 01 + 0.22 
z ‘ Oi - 0.20 

60.1/65 

9 c+0. 54 
z *°-0.59 

61.6/65 

0.41 

2.98 

NGC4631 XMM4 d 

7.8 

9.50 

261.5/74 

2.9 

0.17 

207.8/74 

- 

- 

NGC4631 XMM5 f 

1.3 

1.03 

641.8/153 

1.3 

199 

659/153 

- 

- 

NGC5204 XMM2 

0.89i“; 4 5 ® 

1 QQ+ 0 - 25 

i.yQ — 0 20 

1 - Dd -0.17 

2- 67l°' 2 ! 

42.37/42 

0-23i!;} 2 

4-05li;|! 

42.2/42 

0.15 

0.41 



41.4/47 

0 49 + 0 - 46 
u - 4Z -0.38 

7.82l?;°? 

1 63+ 032 
J..DO_o 25 

39.4/47 

0.25 

0.69 

M51 XMMl 

i 1+0.30 
li -0.27 

110.5/82 

o-istS:}! 

140.4/82 

0.34 

2.8 

M51 XMM3 

0 «+0-30 
U ‘°-0.40 
0 4 + 0-20 
u *^-0.30 

i- 86 i°:°| 

63.2/72 

0.05t°;° 2 

c 99 + 2.26 
^ ZZ -i.41 

69.2/72 

0.18 

1.1 

M51 XMM4 

l-55l°°| 

34.8/37 

o-oitg;}5 

111 +0.32 
ii - i -0.25 

34.8/37 

0.16 

0.99 


a total column density in units of 10 21 cm -2 

b unabsorbed flux in the 0.3-10 keV band in units of 10” 12 erg cm -2 s~ A 

c unabsorbed luminosity in the 0.3-10 keV band, using the distances quoted in Table 1, in units of 10 39 ergs 
d see appendix; super-soft X-ray source best fit by single-component blackbody 
e absorbtion column density fixed to the galactic column density found in Table 1 
{ source is best fit by a combined power law and vapec model; see appendix 
g question mark denotes an error in fitting the parameter 
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Table 3. XMM- Newton best fit: two-component blackbody and power law spectral fits 


Source 

n H* 

kT (keV) 

r 

X 2 /dof 

0 

* 

10 

cr 

Fx c 

Lx d 

NGC247 XMM1 

4 1+ 1 * 9 

0 . 12 l°;°f 

0-80l“;J| 

1 14+0-07 
A ' A 4 - 0.10 

418 ±2:I1 

86.5/93 

25.7 

6.2 

7.1 

NGC253 XMM 1 

2 7+0.4 
Z ‘ —0.4 

i-74lg;}I 

225.9/230 

36.7 

2.7 

4.5 


7 o+i.i 
*•*-0.9 

2 54+ 027 
Z * O4 -0.22 
0 14+0-05 

^^-o.os 

2 47+ 2 -99 
z - 4 -0.41 

567/580 

44.6 

3.4 

5.7 

NGC253 XMM 2 (obs 2 ) 

2 . 0 ±g| 

n 71+0.10 
U -' A -0.10 
0 7^+0- 13 
u - ' -o.io 
0.6T±O; 4 | 

0-lliao3 

n nq + 0 - 02 
u,uy -o.oi 

460.3/498 

47.1 

1.6 

2.7 

NGC253 XMM3 

3-iio:! 

68.5/82 

23.4 

0.60 

1.0 


o O+0.7 
** Z -0.5 

2 n7 +014 
z * u ' -0.20 

347.4/407 

34.4 

0.80 

1.3 

NGC253 XMM4 

9PJ+10.8 

ZU -7.6 

951 +0.49 
Z ' O1 -0.30 

66.7/57 

6.9 

15 

25 


4-5±J;| 

9 OO + 0.27 

^.oo _ 0 22 
2.43±f;g| 

309.3/291 

12.1 

1.4 

2.2 

NGC253 XMM5 

1 £ S “ f_7 - 2 

0 Qfi+°- 24 
U * y °-0.32 

26.5/23 

5.3 

0.26 

0.43 


4 6 4 " 1 * 1 

4 - D -0.7 

0 16 +002 
u 'io _ 0 03 

1-9+2: it 

223.7/296 

60.1 

1.4 

2.2 

NGC253 XMM 6 

6.3±f-| 

0 . 12 i °;“ 2 

2.26t°;}f 

417.9/407 

17.1 

1.9 

3.1 

NGC253 XMM7 

6 -3±?;f 

0.69l“;ll 

2.40l°;ii 

335.8/339 

21.2 

1.1 

1.8 

NGC300 XMM 1 

i 7+0.20 

A ' — 0.30 

0.9S±g;lJ 

0 a-i +0.06 
°-* i -0.26 

443.7/420 

26.1 

1.3 

1.0 

NGC300 XMM2 

3-8111 

0.O9t o o ;°l 

2-87i°i 4 

102.6/97 

31.34- 

1.1 

0.86 

NGC300 XMM3 


0.04l||;“ 

1-98+-; 

87.7/79 

14.2 

1.2 

0.93 

NGC300 XMM 6 

9 o+2.6 

Z *-1.3 

0.84lS;fS 

A 1 O+0.03 
y ‘ A * — 0.02 

4 9+ 1 - 97 

4 - y -0.7 

34.6/35 

13 

0.27 

0.20 

NGC1313 XMM 1 

3-0 ±li 

1 75+0.14 
A -*°- 0.11 
9 97 + 0 . 1 0 

z - z/ -0.14 

194.1/201 

35.4 

0.64 

1.3 

NGC1313 XMM2 

3-liol 

n i fi+ 0 04 
U - iO -0.02 

425.2/419 

38.9 

2.0 

4.2 

NGC1313 XMM3 

£ 9+0.8 
D ‘ Z -0.6 

o-ulHl 

2 76+0. 10 
Z * /D -0.11 

441.7/424 

336.6 

10 

22 

IC0342 XMM3 

q 7+1.8 

0-091H? 

2 6Q+0- 16 
Z-D - 0.23 

129.5/107 

56.3 

31 

56.4 

NGC1705 XMMl 

H 9Q+0-39 
U>zy -0.27 

1 01 +CU1 
A ‘ UA -0.29 

0 23 +01 ° 
u * z *-0.ll 

9 0-1 +0.89 
^•^1-0.48 

53/85 

8.9 

0.10 

0.44 

NGC1705 XMM2 

0 Qft +°* 97 
O.yo_ 0 32 

1 60 +1 ' 97 
i.OU_o,27 

85.5/74 

6.5 

0.09 

0.27 

NGC1705 XMM3 

0.931?' 51 * 

1 n 7+0.20 
A * U '-0.15 

9 90+O.7O 
ZZ *-0.56 

69.8/65 

11.1 

0.15 

0.48 

MRK71 XMMl 

i 4+°- 9 
i - 4 -0.6 

2.3i;i 

i q+ 0.8 
i - 8 -0.6 

1 45+1-94 
i - 4O -0.38 

9 57+1-16 
Z0{ -0.90 

52.6/54 

4.1 

0.24 

0.33 

NGC2403 XMMl 

0.661qJ! 

9 -j O+0.41 
^■ li5 -0.59 

l-95io.il 

81.4/79 

10.8 

1.99 

3.1 

NGC2403 XMM2 

O- 62 T 0 if 

0.74l°; 2 i 

163.1/151 

16.4 

1.0 

1.6 

NGC2403 XMM3 

1 7 + 11 
A * — 0.8 

2 -i3il11 

84.2/105 

8.4 

0.64 

1.1 

Holm II XMMl (obs 1) 

A,D - 0.2 

0 4+°- 5 
U - 4 -0.3 

9 oc + 0.05 
Z '*°-0 .11 

997.5/976 

136.7 

12 

10 

Holm I XMMl 

1 07 + 0.66 

1 -0.89 

2 46+0-44 

z - 1D -0.40 

2 70+° 04 
Z *' U -0.05 

97.4/93 

5.4 

0.6 

0.93 

M81 XMMl 

o /?+0.07 
**°-0.07 

O-Ollaw 

1316.8/1243 

533.1 

5.0 

7.8 


3 *>+ 0A 
7 4+0.5 

-0-T 

1 1 q+0-13 
A * A *-0.14 

2-34±g;|| 

203.5/204 

21.4 

4.8 

7.4 

M81 XMM2 

a -■ +0.004 

O.l_ 0 .004 

2.87il;il 

833.9/616 

524.3 

13 

22 

M81 XMM3 

3 t+ 2 - 4 
**'-2.1 

o.nifoi 

1 fiQ+0-27 

it>y -0.33 

77.1/78 

4.25 

0.53 

0.82 

M81 XMM4 

1 1 +1 - 6 

0.62±° if 
o-i+fos 

2 31 + 1 ' 22 

Z -* A -1.05 

48.9/50 

28.2 

0.43 

0.70 

M81 XMM5 

n 1 c+0.69 
U - A ^-0.13 

1 2fi +0-22 
i - ZO -0.20 

89/80 

8.5 

0.38 

0.59 

Holm IX XMMl 

9 -I +0.2 
Zi -0.2 

1 79+O.O4 
1 • * z — 0.03 

866.6/878 

134.3 

10 

16 

NGC4214 XMM2 

i-b±S:S 

osiiglf 

0 Q5+1-81 

** yo - 1.05 

46.4/44 

4.5 

0.4 

0.35 

NGC4258 XMMl 

0.38+°| 6 

0 54+°- 17 
u, ° -0.08 

1 M+°' 4 
A,OA -0.4 

91.1/76 

10.3 

0.34 

2.1 

NGC4258 XMM 2 (obs 1) 

i q+ 2 - 4 ' 

i<y -0.4 

0-78±g il 

2 . 02 if | 5 

73.4/61 

24.1 

0.31 

1.9 

NGC4395 XMMl 

9 a+0.08 

Z ' U -0.07 

0 5 +0 - 9 
u -°-0.3 

Q- 14 io:o 2 

l-ioloiw 

n 1 r +0.03 
U.lo_ 0 03 

n 77+O-O8 
U *"-0.095 

3 44+0-54 
** 44 -0.56 

168.2/154 

26.9 

1.4 

2.7 

NGC4395 XMM3 

9 aa+ 1 .05 
z -^-0.77 

52/56 

3.9 

0.29 

0.56 

NGC4449 XMM3 

o c + 1-3 
*°-0.9 

9 59+0.36 
ZOZ _0.39 

119.9/87 

34.1 

1.1 

1.3 

NGC4490 XMMl 

c o+2-96 
°-°-2.96 

2 SQ+ 1 77 

z -oy_o. 6 i 

66.5/63 

35 

0.88 

6.4 

NGC4490 XMM 2 

4 4+ 1 - 9 
4 - 4 -1.9 

0.60±°ff 

9 -i o+0. 50 
ZAO -0. 70 

42.4/54 

7.1 

0.65 

4.7 

NGC4490 XMM3 

■i 0+9.6 
A *-2.5 

0.09l“;“ 

0 9 i +0.52 
* ZA -0.17 

72.1/78 

4.6 

12 

87.4 

NGC4631 XMMl 

o+0.9 

*-0.5 

0 - 1 +aol 

9 1 9+0.03 
ZiZ -0.02 

371.3/345 

12.1 

0.96 

6.5 


Table 3 — Continued 


Source 

n // a 

kT (keV) 

r 

X 2 /dof 

Ax 2b 

Fx c 

Lx d 

NGC4631 XMM2 

o o+i- 4 
0.3 

o.i8tS;S 

i-8otg:g 

107.4/97 

12.1 

0.25 

1.7 

NGC4631 XMM3 

i-itoi 

1-01 toil 2 

2-45iJ. 62 

127.1/96 

18.9 

0.15 

1.0 

NGC4945 XMM1 

3.5±f;J 

0 77+ 0 - 27 
u -' '-0.10 

1 - 60 to.31 

96.1/120 

20 

0.59 

0.68 

NGC4945 XMM2 

o o+l- 1 
6Z -0 .7 

1 1^4-0.28 
1 - iO -0.33 

i-80±S:lo 

105.8/113 

8.7 

0.66 

0.76 

NGC4945 XMM4 

4 0 +2 ‘° 

0 10 
U-Oi_o.io 

o qo+1-06 

z ‘ 8z -0.58 

58.4/60 

17.1 

0.38 

0.44 

NGC5204 XMM1 

1 5+ 02 
i,o -0.2 

i - zD -0.09 

1 36 +0 ' 2 
l-'SO— o 1? 

2 qfi+ 0 - 07 
z ' yo — 0.08 

551.1/559 

70.3 

3.0 

8.3 


2 l+°- 4 
zi -0.2 

q ae+0.14 
o.zo_q 12 

566.8/496 

60.1 

5.3 

15.0 

MSI XMM2 

i O+0.6 

— 0.5 

0.26lJ;g 

1 qa+o- 61 
4 -ou_ 0 92 

70.7/68 

4.5 

0.36 

3.0 

MSI XMM5 

10.41J;? 

0.078l°;°i 

2 26 +0 ' 26 
zz »- 0.25 

59.8/70 

196.2 

220 

1900 

M51 XMM7 

2.8±i;i 

o.iotS;g 

1 Q7+O 48 
i ‘ y ' -0.30 

31.7/29 

6.1 

0.26 

1.6 

M83 XMM1 

i 4+0.8 

0.50i°;^ 

A AA + 0.31 
z - zz -0.30 

206.1/209 

4.9 

0.53 

2.5 

M83 XMM4 

0.52j;° 16 

a c^-i“0-06 

0.563o.o7 

i-52t§:S 

95.1/89 

12.4 

0.2 

0.92 

M101 XMM1 

0 22 +9 ' 12 
U ’ zz -0.15 

0 21*h°‘ 03 

U ' zi -0.04 

1-42^0.05 

249.9/231 

53.1 

0.45 

2.9 

M101 XMM2 

1 fi+0-46 
i * c> -0.21 

0 76+ 014 
u - ' D _0.10 

l-SStoIii 

251.6/261 

37.2 

0.7 

4.6 

M101 XMM4 

1 «+ 0 - 17 
1 0.15 

0 54+°- 11 
u *°^-0.07 

o o<5 + 0.12 
z ‘ zz -0.08 

158.2/138 

7.5 

0.34 

2.2 

NGC5408 XMM1 

o.9±8:5S 

0 14 +0 - 01 
U * i4 -0.01 

a 71 4-0.16 
' A -0.20 

316.4/337 

80.4 

3.97 

10.9 

CIRCINUS XMM1 

l0.i±i;| 

D in"*” 0 01 
U ‘ iU _0 .01 

a oa+0 08 
z -* 3U -0.08 

749.4/861 

13.5 

12 

23 

CIRCINUS XMM2 


n ro+0 03 
u - oo -0.03 

4 -71^049 

438.5/430 

79.4 

5.6 

10.7 

CIRCINUS XMM3 

13.5 ill 

o fi7+0.10 
U '°' -0.08 

5-77±l;! 4 

269.3/260 

15.9 

7.6 

14.5 


a total column density in units of 10 21 cm" 2 

improvement in x 2 over the single-component power law model 

c unabsorbed flux in the 0.3-10 keV band in units of 10“ 12 erg cm" 2 s -1 

d unabsorbed luminosity in the 0.3-10 keV band, using the distances quoted in Table 1, in units 
of 10 39 ergs -1 

Question mark denotes an error in fitting the parameter 
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Table 4. XMM-Newton best fit single/ two-component spectral fits 


power law 


blackbody and power law 


Source 

n tf a 

r 

X 2 / dof 

nj/ a 

kT (keV) 

r 

X 2 /dof 

Fx c 

Lx d 

NGC300 XMM5 

?+? e 
• _? 

9 n+0.15 

z - i -o.n 

49.3/55 

0 4i+° 60 
u - 4i ~0.30 

1 06 + 0 - 37 

i.uo_o 20 

9 * 70 + 0.61 
Z> '°-0.65 

47.3/53 

0.17 

0.13 

Sextans A XMM 1 

0 . 181 °;!* 

9 9 C+ 0.12 

Z - ZO -0.07 

271.4/275 

0 4 +0,7 

l-05l^ 7 

9 c+ 0-8 
z -°- 0.2 

269.1/271 

0.60 

0.14 

IC2574 XMM3 

0.i5l°;“ 

2 43 + 0 - 27 

40.3/49 

1 1 + 0.6 
i - i - 0.8 

0.87ljj° 

^•OO_ 0 go 

38.4/47 

0.36 

0.56 

NGC4736 XMM 1 

0-95±g;| 

9 pk9+0.26 

Z - UZ -0.25 

62.8/53 

6-311;? 

0 - 081 °;°* 

0 -lllo.Ss 

9 41 + 0 - 34 
Z ' 41 -0.27 

54.9/51 

8.1 

17.9 

NGC4945 XMM3 

o o + l-3 

i 09 +O.I 2 
iSZ - 0.20 

30.8/30 

6-8 It! 

2 03 +O ' 22 

Z ' UO -0.69 

29.7/28 

0.95 

1.09 

M51 XMM 6 

2-0 t° 0 f 2 

2 ^o -1-0 ' 33 

Z ’ OU -0.25 

40.97/43 

8.2 tii 
1 6 +0 ' 3 

0.08l°;°| 

0 n + 0 - 37 
°* u -0.43 

36.9/41 

5.6 

35 

M 101 XMM3 

1 5 + 04 
lo -0.3 

2 70+°- 21 
Z *' U -0.17 

130.3/133 

0 64 +0 ' 14 
U - D 4 - 0.12 

2.93l°;* 

129.9/131 

0.56 

3.7 

M101 XMM5 

1 3 +0 ' 2 
— 0.3 

9 90 + 0.12 
Z ‘ Z 8 - 0.11 

47.9/46 

1 3 + 1,2 
0.2 

0 1 fi+ 005 
U iy -0.06 

1.951q *2 

45.1/44 

0.13 

0.85 


a total column density in units of 10 21 cm -2 

improvement in \ 2 over the single-component power law model 

c unabsorbed flux in the 0.3-10 keV band in units of 10" 12 erg cm -2 s _1 for combined fit 

d unabsorbed luminosity in the 0.3-10 keV band, using the distances quoted in Table 1, in units of 10 39 erg s -1 

s question mark denotes an error in fitting the parameter 
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Table 5. Bolometric Luminosities of ULX sources 


Source 

r a 

Supper 

Lbol b 

M E dd c 

NGC247 XMM1 

13.4258 ' 

7.07734 

54.4411 

NGC253 XMM1 

9.31469 

2.44574 

18.8134 

NGC253 XMM2 

4.3701 

2.15292 

16.561 

NGC253 XMM6 

5.05828 

3.92514 

30.1934 

NGC1313 XMM3 

37.0364 

27.9692 

215U48 

NGC1313 XMM4 

1.50345 


115.65 

IC0342 XMM1 

14.1215 


1086.27 

IC0342 XMM2 

19.8129 


1524.07 

IC0342 XMM3 

114.015 

95.4068 

733.899 

NGC2403 XMM1 

4.1497 

2.14873 

16.528 

NGC2403 XMM4 

0.57068 


43.898 

Holmberg II XMM1 

0.88906 


68.3893 


16.8335 

11.4543 

88.1103 

Holmberg I XMM2 

10.5158 


808.908 

M81 XMM1 

15.7004 

3.17932 

24.4563 

Holmberg IX XMM1 

31.0582 

28.1445 

216.496 

NGC4214 XMM1 

0.26699 


20.5379 

NGC4258 XMM3 

0.46503 


35.7715 

NGC4395 XMM1 

9.04609 

2.94683 

22.6679 

NGC4449 XMM1 

2.11312 


162.547 

NGC4449 XMM2 

2.48586 


191.22 

NGC4490 XMM1 

16.8513 

3.21972 

24.7671 

NGC4490 XMM2 

7.36612 

4.51554 

34.7349 

NGC4490 XMM3 

240.653 

176.04 

1354.15 

NGC4490 XMM4 

1.66829 


128.33 

NGC4490 XMM5 

12.7136 


977.971 

NGC4631 XMM1 

10.6527 

8.59661 

66.1278 

NGC4736 XMM1 

31.6561 

27.3664 

210.511 

NGC5204 XMM1 

22.4756 

2.20492 

16.9609 

NGC5204 XMM2 

5.57769 


429.053 

M51 XMM1 

4.76208 


366.314 

M51 XMM2 

4.25898 

3.57502 

27.5001 

Mol XMM3 

2.10133 


161.641 

M51 XMM4 

2.56064 


196.972 

M51 XMM6 

46.7642 

39.5189 

303.991 

M101 XMM1 

8.04916 

7.68224 

59.0942 

M101 XMM2 

7.54268 

4.96709 

38.2084 

M101 XMM3 

10.9792 

1.03659 

7.97381 

NGC5408 XMM1 

20.9211 

11.5369 

88.7455 

Circinus XMM1 

70.3579 

56.1033 

431.564 

Circinus XMM2 

208.746 

0.69929 

5.37916 

Circinus XMM3 

771.157 

0.212957 

1.63813 


a upper limit on the bolometric luminosity, determined 
with an exponential cut-off in the power law at high energy 
(see text) 

b boIometric luminosity estimate for high-state ULXs 
where the power law is cut at twice kT{ n (see text); units 
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for both luminosity measurements in 10 39 erg 

c mass computed for objects radiating at 0.1 xL^dd (low- ‘ 
state objects) or L^dd (high-state objects; using L ^/) f in 
units of M© 



Table 6. Blight Point Sources examined 
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NGC4395 XMM2 12 25 25.3 +33 36 46.4 392, 252, 208 2.6, 1.5, 1.2 - sa? 

NGC4395 XMM3 12 25 32.6 +33 25 27.9 763, 277, 278 5.2, 1.5, 1.6 - out? 

NGC4395 XMM4 12 25 42.7 +33 40 00.1 516, 60, 40 3.1, 0.3, 0.2 - out 

NGC4449 XMM1 12 28 18 +44 06 30.9 1409, 608, 593 10.6, 3.6, 3.5 XRB?; NGC4449 [RW2000] X-7; Source 27 13 sa 
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three other sources near the nucleus are unresolveable in the XMM obs., but seen by Chandra (Eracleous et al. 2002) 
Unable to model sprectrum due to an error in x 2 * s Pace 

bright sources are too close, but resolvable by Chandra (Summers et al. 2004) 
alternate ID from Jenkins et al. (2004) 
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Table 7. Bright, Identifiable Background and Foreground Sources 


Galaxy RA (h m s) Dec (o i n) Identification 


NGC 247 

0 46 51.7 

-20 43 30 

NGC 300 

0 55 26.7 

-37 31 25.6 

NGC 625 

01 

34 

42.4 

-41 36 15.2 

NGC 1569 a 

04 

31 

16.9 

+64 49 50 

NGC 1569 

04 

31 

14.2 

+64 51 07.9 

NGC 1569 

04 

31 

25.4 

+64 51 53.8 

NGC 1705 

04 

54 

01.2 

-53 21 12.3 

NGC 2403 

07 

35 

09 

+65 40 27.5 

NGC 4258 

12 

18 

08.9 

+47 16 08.3 

M83 

13 

36 

45.6 

-29 59 13.9 

M83 

13 

36 

13.9 

-29 56 13 

NGC 5253 

13 

39 

50.6 

-31 34 11.1 

M101 

14 

02 

30 

+54 21 18.2 

NGC 5408 

14 

03 

27.5 

-41 25 18.5 


QSO B044-2059 
HD 5403 (Star) 

QSO B0132-4151 

CXOU J0431 16. 8+644950 (Star) 

CXOU 043114.0+645107 (Star) 

CXOU 043125.1+645154 (AGN) 

WGA J0454. 0-5320 (M star or elliptical galaxy) 
HD 59581 (Star) 

QSO J1218+472 

2MASX J 13364579-2959122 (Galaxy) 

RX J133615-2957.8 (Galaxy) 

CD-30 10790 (Star) 

[WIP99] H13 (Star) b 
(Star) 


identification for objects in NGC 1569 from Martin, Kobulnicky, & Heckman (2002) 
b confirmed by K. Kuntz using HST ACS 
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Table 8. Spectral Simulations 


ID 

exposure time (s) 

power law x 2 

bremsstrahlung x 2 

blackbody and power law x 2 


10% exposure time 

NGC 247 XMM1 

1130 

96.4/94 

223,75/94 

72.7/91 

NGC 253 XMM1 

10778 

655.2/581 

667.8/581 

648.7/578 

NGC 1313 XMM1 

2789 

188.7/201 

189.08/201 

187.2/197 

NGC 1313 XMM2 

2788 

425.7/419 

431.7/419 

423.5/415 

IC0342 XMMl 

511 

164.4/183 

208.7/183 

126/179 

IC0342 XMM2 

511 

236.4/107 

189.9/107 

160.9/103 

Holmberg II XMM1 

982 

1034.6/976 

1152.1/976 

971.2/972 

M81 XMM1 

9082 

1314.4/1244 

2010.9/1244 

1246.1/1241 

Holmberg IX XMM1 

709 

898.7/882 

910.7/882 

881.8/878 

Circinus XMM2 

9005 

395.3/431 

388.7/431 

385.6/428 

5% exposure time 

NGC 247 XMM1 

565 

107.9/94 

146.4/94 

77/91 

NGC 253 XMMl 

5389 

516.3/581 

519.2/581 

516.3/578 

NGC 1313 XMMl 

1394 

125.8/201 

125.9/201 

122.9/197 

NGC 1313 XMM2 

1394 

273.3/419 

273.3/419 

263.3/415 

IC0342 XMMl 

256 

141.1/183 

157.7/183 

125.5/179 

IC0342 XMM2 

256 

122.6/107 

99.8/107 

148/103 

Holmberg II XMMl 

491 

908.3/976 

981.8/976 

871.3/972 

M81 XMMl 

4541 

1086.8/1244 

1426.2/1244 

1060.7/1241 

Holmberg IX XMMl 

354 

663.6/882 

635.9/882 

629.2/878 

Circinus XMM2 

4502 

258.7/431 

256.9/431 

255.1/428 
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Table 9. XMM-Newton power law fit for best fit two-component spectra 


Source 

n 

r 

X 2 /dof 

^x c 

Lx a 

NGC247 XMM 1 
NGC253 XMM 1 

9-5111 

3.4±g;| 

8 52+ 1 - 14 
o-O ^_ 0 91 

1 77 +O.O 6 
A *' *-0.06 

112 . 2/95 

262.6/232 

1900 

3.1 

2200 

3.6 


Q + 0 - 4 

°- y -0.4 

1 qq+0.05 
05 

611.6/582 

2.9 

3.3 

NGC253 XMM2 (obs 2) 

O q+ 0.1 
Z * Z - 0.1 

2 o ^"*" 9 ' 04 

507.4/500 

1.2 

1.4 

NGC253 XMM3 

q q+0.5 
°- y -0.4 
4 0 +0 ‘ 4 

tim 

9 17 + 0.14 
Z ‘ A '- 0.11 

91.9/83 

0.73 

1.2 


0 Ofi - * -0 - 09 
^•vO_o.07 

2.09±g;|| 

2.091“;}® 

381.8/409 

0.98 

1.6 

NGC253 XMM4 

73.6/59 

0.52 

0.85 


l- 2 ±g;| 

321.4/293 

0.29 

0.48 

NGC253 XMM5 

1 7 + 1 - 2 

3 4 + 0 I 

1 c S4 +°- 22 

31.8/25 

0.32 

0.53 


9 17+0.7 
Z>i '-0.7 

283.8/298 

1.1 

1.3 

NGC253 XMM 6 
NGC253 XMM7 

0 q+0.3 
7 1 +0.7 

9 91 +0-84 
Z -^ i -0.80 
9 1 r+ 0.11 
z - 10 - 0.11 

435/409 

357/342 

0.93 

1.3 

1.5 

2.2 

NGC300 XMM1 

0 Q 7 +OII 
U * y '- 0.11 

9 £'7 + 0.06 
Z * D ‘ _0.06 

469.8/422 

0.81 

0.6 

NGC300 XMM2 

1 7+0.40 

q 9q+0.32 
J ’ ZU -0.24 

133.98/99 

0.27 

0.21 

NGC300 XMM3 

0 4 +O .8 
5 , > 0.6 

i- 86 ±g;i| 

101.9/81 

0.19 

0.15 

NGC300 XMM 6 

?+? 
' — ? 

2.05i«;” 

47.6/37 

0.06 

0.15 

NGC1313 XMM 1 

1 tr+ 0-2 

A ’°- 0.2 
O q+0.16 
0.16 

1 01 +0.08 
l-ol_0.09 

219.8/203 

0.42 

0.88 

NGC1313 XMM 2 

2 4 .Q + 007 
Z ‘ 45 — 0.06 

464.1/421 

1.9 

4.0 

NGC13X3 XMM3 

3-«±g:i 

q 9+0.09 
^-O.OQ 

778.3/426 

3.3 

6.9 

IC0342 XMM3 

3-s±8:S 

9 CQ + 0.15 
Z '^ 5 -0.14 

185.8/109 

1.7 

3.1 

NCC1705 XMM 1 

0.3 e 

1 93+°- 11 
i - J, 5 -o.io 

9 i 9 + 0.26 
Z * iZ -0.15 

61.9/88 

0.12 

0.37 

NGC1705 XMM2 

1 4 + 0.45 
A * 4 -0.41 

91/76 

0.078 

0.24 

NGC1705 XMM3 

n fi+0-36 
U - D -0.40 

1 36 +012 
r.OO__ 0 .i3 

80.9/67 

0.17 

0.53 

NGC2403 XMM1 

0 q+0.61 

^-o.ss 

9 1 k+0-16 

Z * iO -0.15 
2 07+0- 11 
z,y ' — 0.11 

92.2/81 

2.2 

3.6 

NGC2403 XMM2 

9 7+0.37 
Z -'-0.34 

179.5/151 

1.3 

2.0 

NGC2403 XMM3 

i q+0.40 
i ' y -0.36 
1 r+0.07 
i -°-0.07 

1 Q7+0 - 14 
A * y — 0.13 

92.6/107 

0.81 

1.3 

Holmll XMM1 (obs 1) 

2 fil + 004 
Z ' Di -0.04 

1134.2/976 

12 

10 

Holm I XMM1 


2 04 +0 14 

Z ' U4 -0.07 

102.8/95 

0.48 

1.7 

M81 XMM 1 

0 9+0.07 
^■ Z — 0.07 

2 OQ ^ 0 02 
Z ' Uy — 0.02 

1849.9/1245 

4.5 

7.0 


0 q+ 0.3 
^* u -0.3 

1 7Q+0.07 
1 * /y -0.07 

224.9/208 

4.3 

6.7 

M81 XMM2 

7.3 

6.13 

1358.2/618 

48.5 

75.2 

M81 XMM3 

0.97l°1f 

x -58tg;i« 

81.35/80 

2.5 

3.9 

M81 XMM4 


0 . 88 lg;U 

66.4/52 

0.35 

0.54 

M81 XMM5 

l-Olol 

1 7 +O.O 8 
A * —0.08 

1 co+ 0.11 

L ‘ bZ -o.n 

1 84 +0 * 03 
i * 54 -0.03 

97.5/82 

0.44 

0.68 

Holm IX XMM1 

1000.9/882 

9.4 

15 

NGC4214 XMM2 

n 2+ 0 - 5 
U * z - 0.2 

2 03 +0,43 
z - U6 ~0.2S 

50.9/46 

0.16 

0.14 

NGC4258 XMM1 

i-«±8:5 

1 q+0.14 

i * y -0.13 

101.4/78 

0.06 

0.04 

NGC4258 XMM2 (obs 1 ) 

3.5l°;® 

l.sslgH 

97.5/63 

0.43 

2.6 

NGC4395 XMM 1 

q 7 + 0-5 

4-93ig;i^ 

195.1/156 

7.2 

14 

NGC4395 XMM3 


i-oO _ 0 09 

55.9/58 

0.25 

0.48 

NGC4449 XMM3 

q q+0.5 

q q£+0.29 

<5"5 d _0.23 

9 tcq+O - 17 
^•9O _ 0 16 

154/89 

1.3 

1.5 

NGC4490 XMM 1 

0-83±g;l| 

101.5/65 

1.2 

8.7 

NGC4490 XMM2 

fi Q + l-3 

2 qfi + 019 
l ' 6b -0.\7 
9 Qc;+ 0- 24 

z.yo _ 0 .20 

49.5/56 

0.92 

6.7 

NGC4490 XMM3 

Q 4 + 1 ' 5 
y ‘ 4 - 1.2 

76.7/80 

1.5 

11 

NGC4631 XMM1 

9 q+0.16 

Z ' 6 -0.1S 

9 1 q+0.06 

Z * 1,5 -0.05 

2 Ol +014 
Z ' Ui - 0.12 

383.4/347 

0.76 

5.1 

NGC4631 XMM2 

1 Q+ 0 - 4 

1 —0.3 

119.5/99 

0.23 

1.5 
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Table 9 — Continued 


Source 


r 

X 2 /dof 

Fx c 

Lx* 

NGC4631 XMM3 

0.63±g'| 

i.oo_q 08 

l*oo_0.lO 

A - & °-o.io 

9 ^Q^" 9,19 
z.oy_ 0 I? 

9 11 +0-04 
Z ' 1 — 0.04 

146/98 

0.15 

1.0 

NGC4945 XMM1 

3 4+ 0 - 6 
d -^-0.5 

116/122 

0.9 

1.0 

NGC4945 XMM2 

114.5/115 

0.71 

0.82 

NGC4945 XMM4 

5-2to'.r 

75.5/62 

0.49 

0.56 

NGC5204 XMM1 

0.6lt°;l 

621.4/561 

2.0 

5.5 


1 1+ 0 * 1 

O 41 + 0 07 

2 ^0+°‘ 22 
Z .OU_ 0 20 

626.9/498 

3.0 

8.3 

M51 XMM2 

9 q+0-50 

^-0.30 

75.2/70 

0.52 

3.3 

M51 XMM5 

2.7 

3.08 

256.0/72 

0.43 

2.7 

M51 XMM7 

rt tr+0.39 
u *°-0.46 

1.95±g;f| 

37.8/31 

0.11 

0.66 

M83 XMMl 

i q+0.34 
i - y -0.31 

4.4±J;| 

9 09+0.13 
Z ” 5Z — 0. 12 

210.9/211 

0.64 

2.8 

M83 XMM4 

2-30±g;?g 

1 QO+0 -08 

i.yo_ 0 os 
1 as”*" 9 ' 97 

i-oO_ 0 _07 

107.5/91 

0.4 

1.8 

M101 XMMl 

0.56i°;i= 

303/233 

0.45 

2.9 

M101 XMM2 

9 9+0.25 
Z * Z -0.23 

288.8/263 

0.81 

5.3 

M101 XMM4 

9 9 + 0.45 
Z ‘ Z — 0.42 

9 O^+O ^Q 
z * zo -0.17 

165.7/140 

0.38 

2.5 

NGC5408 XMMl 

1 fi - * -9 ' 2 

o c 7+0. 12 
6 - 0{ -0.11 

396.8/339 

7.04 

19.4 

CIRCINUS XMMl 

7-6±g;| 

9 i r+0.05 
z ' iO -0.04 

762.9/863 

4.6 

8.8 

CIRCINUS XMM2 

ll-7tg:r 

Q 40+°- 13 

9 ry+0.40 
Z -°'-0.17 

517.9/432 

2.7 

5.2 

CIRCINUS XMM3 

9-ot?:o 

285.2/262 

0.32 

0.61 


a total column density in units of 10 21 cm" 2 

b improvement in \ 2 over t ^ ie single-component power law model 

c unabsorbed flux in the 0.3-10 keV band in units of 10 -12 erg cm -2 s" 1 

d unabsorbed luminosity in the 0.3-10 keV band, using the distances 
quoted in Table 1, in units of 10 39 erg s -1 

e absorption level frozen at approximate galactic level 
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Table 10. Best-Fit Absorbed Comptonization Model Parameters 


ID 

nH a 

kT b 

tau c 

X 2 

Fx d 

NGC 253 XMM2 

i q + 0.08 

1 2S -,_0 ' 13 
n QR + 0.16 
U * y °- 0.15 

0 Qfi+o.13 
U ' y0 -0.16 

i,Zi -0.31 

19.59^ 

464/498 

1.47 

NGC 2403 XMM1 

25-4tS:S 

82.8/85 

1.4 

NGC 4490 XMM1 

4-T±l \ 

27.0lJ 8 9 3 

66.5/64 

0.66 

NGC 4490 XMM2 

5-o±i:S 

18.8+1? 

45.7/55 

0.67 

M101 XMM2 

1 e + 0.26 

1>D -0.24 

1 24 +0 * 09 

23.3±|J 

256/262 

0.65 

M101 XMM3 

1 -i +0.43 
A ‘ — 0.40 

1.13 

15.2 

128/132 

0.41 

Circinus XMM2 

a 0 + 1-3 
bii -0.9 
6 7+ 1 - 6 

0-62±°;°! 

29.7+J 

437.2/430 

0.5 

Circinus XMM3 

0-93±°;l® 

23.1±|^ 5 

273.1/261 

0.17 


a total column density in units of 10 21 cm" 2 
b temperature in keV 
c optical depth 

d unabsorbed flux in the 0.3-10 keV band in units of 10“ 12 erg cm -2 s _I 
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Table 11. XMM-Newton Galaxy Observations 


Galaxy 

S 60 (Jy) 

Sioo (Jy) 

Ffir* 

L FIR b 

No. of ULX 

NGC247 

7.93 

27.32 

0.602 

0.687 

1 

NGC253 

998.73 

1861.67 

55.92 

93.10 

3 

NGC300 

23.08 

74.45 

1.688 

1.324 

0 

NGC625 

5.09 

9.08 

. 0.280 

0.230 

0 

NGC1313 

35.97 

92.00 

2.329 

4.845 

2 

IC0342 

255.96 

661.68 

16.66 

30.32 

3 

NGC1569 

45.41 

47.29 

2.072 

0.635 

0 

NGC1705 

0.970 

2.580 

0.064 

0.199 

0 

MRK 71 

3.51 

4.67 

0.173 

0.239 

0 

NGC2403 

51.55 

148.49 

3.547 

5.378 

2 

Holmberg II 

1.15 

2.62 

0.070 

0.061 

1 

Holmberg I 





1 

M81 

44.73 

174.02 

3.647 

5.655 

1 

M82 

1271.32 

1351.09 

58.35 

106.2 

1 

Holmberg IX 





1 

Sextans A 

0.255 

0.674 

0.017 

0.004 

0 

IC 2574 

2.41 

10.62 

0.212 

0.329 

0 

NGC 4214 

17.87 

29.04 

0.947 

0.826 

1 

NGC 4258 

21.60 

78.39 

1.690 

10.48 

1 

NGC4395 

4.21 

12.90 

0.299 

0.573 

1 

NGC4449 

37.00 

58.28 

1.937 

2.199 

2 

NGC4490 

47.79 

85.94 

2.636 

19.19 

5 

NGC4631 

82.90 

208.66 

5.324 

35.83 

1 

NGC4736 

62.41 

135.34 

3.734 

8.261 

4 

NGC4945 

588.11 

1415.5 

36.95 

42.49 

0 

NGC 5204 

2.33 

5.35 

0.143 

0.395 

2 

M51 

108.68 

292.08 

7.213 

44.74 

5 

M83 

266.03 

638.63 

16.69 

76.79 

0 

NGC5253 

30.00 

30.92 

1.365 

16.72 

0 

M101 

88.04 

252.84 

6.048 

39.63 

4 

NGC5408 

2.825 

2.958 

0.129 

0.356 

1 

Circinus 

248.7 

315.85 

12.06 

23.10 

4 


a flux in units of 10 9 ergcm 2 s 1 
b far-infrared luminosity in units of 10 42 ergs' 1 
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l0 g( f x/ f opt) 


Fig. 2 — The distribution of f x /f op t for the brightest possible optical point source within the 
XMM-Newton error circle. We define f x as the unabsorbed X-ray flux in the 0.3 — 10 keV range and 
fopt as the optical flux obtained from the U filter of XMM s OM (as described in text). These ratios 
do not represent the actual f x / fopt of the sources but are an estimate of the maximum possible 
value. A majority of the sources had no optical point source within the X-ray contour and thus 
have ratios of f x /f 0 pt far higher than those indicated in the plot. 




Spectral Index (f) 


Fig. 3. — Distribution of the spectral indices (r) for low-state (shaded) and high-state (open) 
objects. For Galactic low-state objects, typically T « 2.0, similar to our sample, while the high- 
state objects have a steeper T (McClintock & Remillard 2004). 
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Spectral Index (r) 


Fig. 4. — Relationship of luminosity vs. spectral index for low-state (rectangle) and high-state 
(triangle) objects. As expected from observations of Galactic stellar-mass black hole systems (Mc- 
Clintock & Remillard 2004), the classified low-state ULXs in our sample have, on average, lower 
X-ray luminosities than the corresponding high-state ULXs. We plot the mean values for both 
high-state and low-state objects with errorbars indicating the root mean square deviation. The 
outlying objects with spectral indices greater than 3.5 were not included in the mean or deviation 
calculations. 
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41 


40.5 


3 

I 

40 


39.5 

-1 - 0.5 0 

log(kT) 

Fig. 5. — (left)Distribution of the blackbody temperature for high-state objects, (right) Relationship 
of blackbody temperature vs. luminosity (in the 0.3-10 keV band) for high-state objects. We see 
two peaks arise in the distribution, one centered around kT« 0.1 and another at kT« 1. The peak 
with a low disk temperature also corresponds to the highest luminosities, suggesting that these may 
be high-state IMBHs. The sources with higher disk temperature also have lower luminosities. The 
spectra of these sources were also well fit, by an inverse Comptonization model (a model succesfully 
used to fit some of the Galactic black hole X-ray binaries in the very high state). 



F”i i 

, , , , r* 

. « « r ' = 

x 

r 





- 

X 



- 

L 


_ 

! x 

- 


- 

- 


- 

X 

X 

- 


- 

■ X 

■ 

X 

X X - 

- 


- 


X 

X 


X 


X 


X 

X 

X 

X 


X 

X 



X ^ 

' 


X 

“ X 

X X 

X - 


X 

X 

■ , 1 

I 1 > 1 . 

. . 1 , ~ 

-1 

-0.5 

0 


log(kT) 




-42- 



Fig. 6. — Color-Color Diagram plotting soft vs. hard colors, as outlined in Done & Gierlinski 
(2003), for low-state (rectangle) and high-state (triangle) ULXs. A large number of our sources lie 
in the same range of this graph as the black hole sources examined by Done & Gierlinski (2003) 
(near the power law distribution, indicated by the solid line). The dashed line represents the color- 
color plot for a multi-colored disk model with different disk temperatures. The sources approaching 
this line were those well-fit by the Comptonization model. Done & Gierlinski (2003) had no black 
hole sources in this region, but atolls and Z-sources, which were also well-fit by Comptonization 
models. 



5 




Fig. 7. — (left) Relationship of the far-infrared luminosity, as an indicator of star formation rate, vs. 
the number of ULXs for each galaxy. If ULXs are associated with star formation, we naively expect 
that the higher the FIR luminosity the more ULXs the galaxy will host, (right) The distribution of 
average number of ULXs / L fir bin for spirals (solid line) follows this expectation. The distribution 
of irregulars (dashed line) is not so easily interpreted. The numbers at the top indicate the number 
of spirals/irregulars in each of the luminosity bins. More irregulars would need to be included in 
this survey for meaningful statistics on this group. 
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log(n H ) 

Fig. 8. — Distribution of the hydrogen column densities of ULX sources. The n# values were 
obtained through spectral fits using the wabs model in XSPEC. Galactic column densities towards 
the host galaxy were subtracted from the spectral fit values. A majority of our ULX sources have 
high column densities (> 10 21 cm -2 ), suggesting that some of this absorption originates with the 
local ULX environment. Bins to the left of the dashed line represent sources with column densities 
very close to the Galactic value and thus a simple subtraction is not statistically representative of 
the true value. 




